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1. Introduction 

in . 

Understanding transport properties in open many-particle systems is one of the main 
goals of contemporary nonequilibrium statistical mechanics. The ultimate goal is 
that of finding the statistical measure for stationary out-of-equilibrium conditions. 
In fact, this would allow evaluating the fluctuations of the relevant macroscopic 
observables (such as the currents) and, possibly, deriving the corresponding transport 
equations, without any ad hoc statistical assumption. In view of the many technical 
difficulties that one typically encounters along this path, it is convenient to start 
investigating simple models, like chains of nonlinear oscillators [HE]- A particularly 
interesting system is the Discrete Nonlinear Schrodinger (DNLS) equation O 2] 
that has important applications in many domains of physics. A classical example 
is electronic transport in biomolecules j5], while in optics or acoustics it describes 
the propagation of nonlinear waves in a layered photonic or phononic media [6l [7] . 
With reference to cold atomic gases, the DNLS equation provides an approximate 
semiclassical description of bosons trapped in periodic optical lattices (for a recent 
survey see [5] and references therein). 



This system is rather interesting since the presence of two conserved quantities 
(energy and number of particles) naturally requires arguing about coupled transport, 
in the sense of ordinary linear irreversible thermodynamics. In fact, in spite of the 
many studies of heat conduction in oscillator chains [TJ [S] [TTJ [T^] , coupled transport 
processess have been scarcely investigated [T3l [14j [15j [16j H3- Furthermore, the 
DNLS equation has been mostly studied because of the breather solutions that can 
spontaneously emerge [18j , while the analysis of its nonequilibrium properties is still 
in its early days [9J [10] . 

The statistical analysis of any system of physical interest requires a proper 
modelling of the interaction with an external reservoir. The reservoir is expected 
to exchange energy and particles with the system until a steady state is reached, 
characterized by the expected temperature and chemical potential. In standard 
oscillator chains, this is typically achieved by introducing Langevin thermostats (e.g., 
see [U12])- In models like the DNLS equation, this option is less straightforward, 
because of its peculiar structure and symmetries. In this paper we augment the 
Hamiltonian equation with a suitable nonlinear damping and a stochastic term to 
model the interaction with one or more reservoirs. To our knowledge, this is the first 
such scheme to be proposed in the literature, at least in the present context, although 
one should mention [19] , where the evolution of a DNLS system has been discussed 
in the presence of small nonlinear damping and a multiplicative noise. 

This general Langevin scheme is then used to investigate transport properties in 
two limit cases, low temperatures and large particle densities, where the DNLS model 
reduces to a chain of coupled oscillators with just internal forces. In fact, it turns 
out that norm (number of particles) conservation is mapped onto the conservation of 
(angular) momentum. Such a relationship helps also to understand the reason why 
DNLS chains are characterized by a normal [10 , rather than anomalous transport, 
as in typical translationally invariant nonlinear systems [1 . The corresponding 
oscillators, in fact, involve phase-like variables such as the XY model: in such models, 
phonons are scattered by the occasional phase-jumps scatter the phonons, thereby 
contributing to a normal diffusion [JJ. Moreover, we find that the chemical potential, 
which is associated with norm conservation, turns out to be equivalent to the rotation 
frequency of the single rotors and the corresponding force that must be exerted by 
the external Langevin reservoir for its thermalization is an additional constant torque. 
Finally, the possibility to map the DNLS equation onto a standard chain of coupled 
(phase) oscillators allows deriving a local microscopic definition of the temperature, 
based on their kinetic energy. 

The paper is organized as follows. In section [2] we recall some basic properties 
of the model at finite-temperature. In section [3] we present the Langevin equations 
and discuss the equilibrium setup as well as the case of two external reservoirs at 
different temperatures and chemical potentials. The rest of the paper focuses on two 
limit cases: we first consider the weakly-nonlinear limit, which can be analysed by 
introducing a suitable harmonic approximation (section [4]), and then pass to the large 
particle-density limit (section [5| , where the dynamics can be mapped onto that of a 
chain of coupled (nonlinear) rotors. A numerical test on the corresponding definitions 
of kinetic temperatures is provided in section [6j while section [7] is devoted to a final 
discussion about the achievements and perspectives of the present study. Finally, the 
appendix contains an alternative derivation of the generalized Langevin equations. 



2. The DNLS model at equilibrium 

The isolated one-dimensional DNLS equation, defined (in suitable units) on a lattice 
of N sites, reads 

iz n = -2\z n \ 2 Z n - Z n +1 ~ Z n -\ , (1) 

where the quantity \z n \ 2 can be interpreted as the number of particles, or, equivalently, 
the mass at site n. Eq. ([T]) can be derived from the Hamiltonian 

N 
H= Y^ (\ Z n\ 4 + Z^Zn+1 + Z n Z? l+1 ) . (2) 

n=\ 

The sign of the quartic term is assumed to be positive, i.e. we consider the case 
of repulsive on-site interaction; the sign of the hopping term is irrelevant, due to 
the symmetry associated with the canonical (gauge) transformation z n — >• z n e l7Tn . 
With this choice of the sign of the nonlinear term and assuming periodic boundary 
conditions, the ground state of the isolated chain is the time periodic solution 
z n = y/a e ^ [ wt + ,^,l ] , where the parameter a fixes the value of the mass density, with 
u = 2(o-l) [2Q]. 

The equilibrium statistical mechanics of the continuous nonlinear Schrodinger 
equation was first studied in [21], while that of its discrete, nonintegrable, version 
([U was developed more recently in [20]. Several conclusions of these studies have 
been extended to a wider class of nonlinear lattices, whose dynamics is described by 
DNLS-like equations [22]. 

An important property of the DNLS dynamics is that it conserves two quantities, 
the energy @ and the total mass 

N N 

^| 2 =2^ (p ™ + 9 ™ ) ' (3) 

n— 1 n— 1 

where we have introduced, for later reference, the notations z n = (p n + iq n )/y/2. 
As a consequence, the equilibrium phase-diagram is two-dimensional, as it involves 
the energy density h = H/N and the mass density a = A/N (in a microcanonical 
description of thermodynamic equilibrium) or, equivalently the temperature T and the 
chemical potential /i (in a grand-canonical description of thermodynamic equilibrium). 
The first reconstruction of the equilibrium phase-diagram was reported in [20] for 
the grand-canonical ensemble, with the help of transfer integral techniques. For our 
purposes, it is sufficient to recall that the positive-temperature states are located in 
the region of the (a, h) plane lying between the ground-state line (see Fig. [1| 

h = a 2 ~2a T = 0, n, = 2(a - 1) 

and the infinite-temperature one (corresponding to random phases) 

h = 2a T = oo, /i = — oo . (4) 

Above the T — oo line, one finds the so-called negative-temperature states |20] , 
In this region, the dynamics of the DNLS equation is characterized by long-lived 
localized excitations (discrete breathers) ^23 . Its thermodynamic properties will not 
be considered here: we refer to |24j and to the literature cited therein for an account on 
the role of these peculiar states on the equilibrium properties of the DNLS equation. 
The nonequilibrium studies discussed in this paper all refer to the positive-temperature 
region. 
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Figure 1. Equilibrium phase diagram of the DNLS model in the (a,h) plane. 
The positive temperature region lies between the ground state (solid blue line) 
and the infinite temperature isothermal (dashed red line). The line at constant 
chemical potential /i = 2 has been reconstructed using the Monte-Carlo stochastic 
thermostats 1101 (purple triangles) and the Langevin scheme (solid black line). 



3. Langevin thermostats 

Being interested in exploring the thermodynamic properties of a Hamiltonian system, 
the first problem we have to deal with is that of introducing suitable thermal 
baths. In separable Hamiltonian systems (i.e., those composed of a kinetic and a 
potential energy) it is sufficient to modify the momentum equation, by adding a 
linear dissipation term of the type —jpt accompanied by a white-noise fluctuation, 
whose amplitude determines the value of the temperature. This simple Langevin-like 
scheme does not work for the DNLS. In fact, one can easily check that in the absence 
of fluctuations (i.e. in the zero-temperature limit) this dissipative dynamics indeed 
converges to a fixed point, that does not correspond to the ground state, which, as 
mentioned in the previous section, is a time-periodic solution |20j . 
The problem can be overcome by adopting the following scheme, 

iz n = (1 + il) [-2\z n \ 2 z n - Z n+1 - Zn-l] + il^z n + y/jT^ n (t) , (5) 

where £„ (t) = £' n + i£^ is a complex, Gaussian, white random noise with unit variance. 
In practice, the above equation corresponds to a series of thermostats all operating at 
temperature T and chemical potential /u, and with the same coupling strength 7 to 
the single DNLS degrees of freedom. This is basically a stochastic, discrete, complex 
Ginzburg-Landau equation. As required for a meaningful reservoir, the dissipative 
term vanishes for the ground state evolution, z n — ^/aexp [i{u>t + 7m)]. 

In the Appendix we show that the Langevin scheme in Eq. (JSJ can be interpreted 
(with the help of suitable assumptions and approximations) as the action of an 
ensemble of linear oscillators. Some physical insight can be, however, gained by simply 
rewriting Eq. ([5]), in terms of the p n q n variables, 



{vl 



iDq-n 



Qn-l -7- 



dp n 



V*yT&{t) 



(6) 



dH u 



Qn = (Pi + ql)Pn + Pn+1 + Pn-1 ~ 7"^ + \/ 2 7 T £n(*) ) 

where and 7J M is the effective Hamiltonian H^ = H — fj,A. The two deterministic 
terms, proportional to 7 in the r.h.s. of the above equations, are gradient terms: 
in the absence of the symplectic component, they would drive the system towards a 
time- independent state with a minimal H^. On the other hand, the addition of the 
symplectic forces ensures a (periodic) motion across the microstates characterized by 
the same H^ value. Accordingly, it is reasonable to conclude that this is the proper 
way to including a dissipation scheme in the DNLS case. The reason why H^ is 
considered instead of H is the presence of two conservation laws: the minimum of the 
energy depends on the mass density a that has to be, on the average, preserved and 
this can be achieved by adding the term — \xA into the effective Hamiltonian, where 
(i plays the role of the chemical potential. The additional presence of the fluctuation 
terms completes the definition of the generalized Langevin equation, that describes 
the coupling of the DNLS chain with a proper stochastic reservoir with temperature 
T and chemical potential /1. 

The Langevin equations (JSJ can be straightforwardly generalized to a non- 
equilibrium setting, by assuming that the single reservoirs operate at different 
temperatures/chemical potentials. In the following we will focus on a typical setup in 
which only the first (z±) and the last (zjv) lattice variables interact with the reservoirs. 
This means that, assuming fixed boundary conditions (i.e. zq — 0), the evolution on 
the leftmost site is rule by the equation 



pi= - (pl + qf)qi - 92 - 7 [(pi + ii)pi +P2- hlpi] + VZjTlZ'i (7) 
<?i = (p? + ql)pi +P2-1 [bi + q\)q\ + <?2 - hlQi] + \J^T L %{ , 

where Tf, and /x^ denote the temperature and the chemical potential of the left 
reservoir, respectively. Analogous equations hold for the right reservoir, which acts 
on the site n = N, where the temperature is Tjj, the chemical potential is /x#, and 
the coupling strength again equal to 7. The rest of the chain follows the Hamiltonian 
evolution ((T|). In order to test the correctness of the method, we have frist simulated a 
chain in contact with two reservoirs operating at the same temperature and chemical 
potential. In Fig. [IJ we compare the isochemical line \x = const, determined from 
the Langevin dynamics (0, with the result with a Monte Carlo thermostat, as the 
one described in [10;. The agreement is well within the statistical fluctuations. 
Moreover, in Fig. [2] we show a typical relaxation process towards an equilibrium 
state characterized by the temperature and the chemical potential imposed by the 
reservoir. In fact, the inset in the same figure shows that the asymptotic state is 
spatially homogeneous as revealed by the flat mass and energy density profiles. 

Although the focus of this paper is on low-temperature transport, it is instructive 
to analyze first the limit of large temperatures. In separable Hamiltonian models, 
large temperatures necessarily correspond to large momenta and large fluctuations 
of the displacement variables. In the DNLS case such conditions occur only if also 
the chemical potential \i is sufficiently large. In this limit, in order to avoid that the 
Hamiltonian dynamics is washed out by the fluctuations, it is necessary to assume a 
vanishing dissipation. In the limit T — ¥ 00, /i — > — 00 (stability is ensured only for 
negative values of n) and 7 — > 0, Eq. ([5]) simplifies to 

izi = ~2\zi\ 2 z 1 - z 2 + i~fii L zi + y/iT L £i . (8) 
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Figure 2. Evolution of T (black solid line) and fi (red dashed line) in a chain 
of N = 100 sites in contact with two Langevin heat baths with TV, = Tr = 1, 
Ml = t^R = 2. The initial state corresponds to a = 4 and h = 15. T and /t 
are measured according to I10| and averaged over running windows of 50 time 
units. In the inset: final spatial profiles of mass (blue dots) and energy (green 
dot-dashed line) densities as a function of the rescaled site index x = n/N. 



Upon introducing T = —J/j,l > and a^ = —Tl/^l, this equation can be written as 

(9) 



i-i 



-2\zi\ 2 zi - z% - iVz\ + volT£i 



This is the "naive" DNLS Langevin equation, which thus we see it corresponds to 
the infinite temperature limit, i.e. it simulates only the points located along the line 
defined by Eq. ((4]). In Fig. [3] we show that this holds true also in the nonequilibrium 
setting, when ol is assumed to be different from clr. In the steady state, the energy 
and density profiles fulfill the relation (j4]) at any point of the chain. Moreover, we 
have checked that the distribution of the local mass is Poissonian, as expected in the 
T = oo case [53] (see the inset of Fig. [3J . 

4. The low-temperature limit 

In order to explore the low-temperature limit, it is useful to express the solution of 
the dynamics (Jl} as 

z n = Z(l + C„) exp[i(ut + </>„ + mr)] . (10) 

Notice that the ground state solution corresponds to £„ = and cj> n = with Z = i/a 
and w = 2(Z 2 - 1). 

By inserting flTOf into Eq. ([1]), one finds that the new variables ( n and <f> n obey 
the following dynamical equations 

0n(l + Cn) = 2£?C„ + 2Z 2 (3C^ + C) - (1 + Cn+l) COS(0„ +1 - <j> n ) 

-(l + Cn-l)cOs(<?(>„-0„_i) + 2 (11) 

Cn = (1 + Cn+i) sin(0„ + i - <p n ) - (1 + Cn-i) sin(0„ - n _i) , 




Figure 3. Noncquilibrium profiles of mass density (upper panel) and energy 
density (lower panel) obtained with the infinite temperature Langevin equation 
(0 and parameters a^ = 0.6, an = 0.5, T = 1. Black solid, red dashed and 
blue dot-dashed lines refer respectively to chain lengths N = 25, 50, 100. The 
profiles are almost linear and h{x) = 2a? (x) along the chain, confirming that the 
thermostats act at T = oo. The inset shows that histogram of the local mass 
\z n ? at n = 20 is Poissonian, as expected. 



where B = 1 + 2Z 2 . In the low-temperature limit, the solution of these equations is 
expected to be close to the homogeneous periodic motion of the ground-state solution. 
Thus, we can assume that C, n <S 1 and (<p n — 4>n-i) < 1. As a result, Eq. (fTTj) simplifies 
to 



(fin = 2-B£n — Cn+1 — Cn-1 
Cn = (<?Wl - 20„ + (j) n -i) 

If one now introduces the new variable 

Pn = %BC, n — Cn+1 — Cn-1 

the Eqs. (fTIZj) can be re-written as 



Pn 



Pn 

2(1 



(12) 

(13) 
(14) 



B){<t>r, 



2(j) n + 4> n -l) - <j) n +2 + 20„ - 4> n 



These equations describe the dynamics of a chain of harmonic oscillators with nearest- 
neighbour and next-to-nearest-neighbour. The corresponding Hamiltonian, 



n h = 






1 2 
~ 2 Pr, 



(1+W 



n+1 



-) 2 



Un+2 



(15) 



is, at leading order in pi and (</>«+ 1 —(fin), fully equivalent to that of the original DNLS 
equation. Its quadratic structure is hardly a surprise, since it essentially originates 
from a parabolic approximation around the minimum of the energy. The normal 
modes, i.e. the plane- wave solutions of Eqs. (|14l) . are the discrete analogs of the 
Bogoliubov modes, in the context of the physics of atomic condensates [26] . 



Passing to thermodynamics, a first interesting implication of the Hamiltonian 
structure in the low-temperature limit is that one can naturally introduce a 
microscopic definition of temperature as the average kinetic energy, i.e. 

T h = J((2BCn ' Cn+l + C«-l) 2 ) , (16) 

where the proportionality constant J is the Jacobian determinant of transformation 
P3[) . which must be included to allow for a meaningful comparison with the DNLS 
model. We have measured Th by numerical simulations of the Langevin scheme 
dchncd in Eq. (JSJ and in Fig. 2] we have compared it with the temperature of the 
bath, Tb = Tl = Tr, for different values of the mass density, a. As expected, Th 
approaches Tb for increasing values of a. In fact, the larger is a, the smaller is the 
relative amplitude of the fluctuations with respect to the ground state. 




Figure 4. Comparison of the harmonic temperature Th, Eq. J16H with the 
reservoir temperature, Tg. The dashed line correspond to 7), = Tg. The 
Langevin heat baths are coupled at the boundaries of a DNLS chain with N = 50. 
Th is measured on a subchain of 30 lattice sites to increase the statistics. The 
values of the Jacobian determinant J, which are the product of the eigenvalues of 
a tridiagonal matrix, have been computed analytically in the N — > oo limit and 
are J = 5.95, 3.73, 2.91, 2.42 for curves from top to bottom. 



Altogether, the original DNLS equation in the low-temperature limit maps onto a 
translationally invariant harmonic chain with nearest- and next-to-nearest-neighbour 
interactions (the former being attractive and the latter repulsive). Moreover, the total 
mass-convervation law of the DNLS maps onto the conservation of the total momentum 
P = ^2p n - The harmonic structure of Hamiltonian (|I5j) implies that the dynamics 
can be decomposed into that of non-interacting phonons. As a result, for what 
concerns transport properties, in the low-temperature regime, the heat conductivity 
should exhibit a divergence in the thermodynamic limit, as expected for any harmonic 
(integrable) model (see |27[|Tj). On the other hand, previous non-equilibrium studies 
of the DNLS model have revealed that the heat conductivity is finite 10 at finite 
temperatures. This is clearly a consequence of the nonlinear terms which break the 
integrability of the dynamics. However, it is not obvious how this can happen in a 
translationally invariant model such as (1151) . since it is known that in the absence of 



pinning forces, the heat conductivity diverges in the thermodynamic limit [2[TT]. The 
problem is solved in the next section, where we discuss the tight relationship with the 
one-dimensional XY model. 

5. Large mass-densities 

Another interesting limit of the DNLS model is that of large a values, that is equivalent 
to small-amplitude fluctuations. This limit can be analyzed by applying a different 
perturbative scheme, which allows preserving relevant nonlinear terms, thus yielding 
more meaningful predictions. The starting point is again Eq. (|f f [) . where we assume 
Z ^> 1 and C, n ~ O(l)- The dynamical equations ([l) simplify to 

4> n = 2Z\ n (2 + 3C 2 , + C)/(l + C„) (17) 

4 = (1 + Cn+i) sinO„ +1 - <p n ) - (1 + Cn-l) shi(</>„ - (j) n -i) . 

It is not clear whether these equations may describe a Hamiltonian dynamics (in 
fact, (f> n and £ n are not canonical variables). If, however, we introduce the additional 
assumption £ n <C 1 and rescale ( n to A„ = 4Z 2 ( n , we obtain, to leading order, 

4>n = Xn (18) 

A„ = 4Z 2 [sin(0 n+ i - (j) n ) - sin(0 n - <j> n -ij\ , 

which corresponds to a system of coupled rotors, i.e. a classical version of the XY 
model in one dimension 28, 29] 1301 E] , with Hamiltonian 

H X Y = Y,Y~^2 Az2c08 ^ n+i ~ <t)n) ' (19) 

n n 

where A„ and 4> n act as a couple of conjugated action-angle variables, the former 
playing the role of an angular momentum. This analogy was already noticed in a 
2d setup in |32j . At variance with the previous section, we have not introduced any 
smallness hypothesis for cf> n — 0„_ i; this is the reason why some nonlinear terms have 
been kept. 

The XY Hamiltonian (|19[) can be obtained more rigorously by applying a sequence 
of canonical transformations to the original DNLS model. Let us first introduce the 
transformation to action-angle variables z n = y/a^ exp[i(tut + <fi n + mi)] and re-write 
the DNLS Hamiltonian (J2J in the form 

H = ^ a n - 2^/a n a n+1 cos(</>„+i - <j> n ) - ua n , (20) 

n 

where the last term in the Hamiltonian follows from the time-dependence of the 
transformation. Taking inspiration from the procedure discussed in the previous 
section, one could then set a n = Z 2 (l + Cn) 2 - This transformation is, however, 
not canonical and this is the reason why Eqs. (TT7|) are not symplectic. One can, 
nevertheless, proceed by defining a n = Z 2 +2Z 2 ( n , a (linear) canonical transformation, 
which approximates better and better the previous expression for increasing values of 
Z. By then recalling that to = 2(Z 2 — 1) and expanding in Z, the Hamiltonian (fT9"|) is 
finally obtained, up to an irrelevant constant shift. 

It is interesting to observe that in this large-density limit the invariance under 
rotation of the DNLS transforms into the invariance under a translation of the angles 
<f> n . Accordingly, the conservation of the total mass A transforms into the conservation 
of the total angular momentum L = J2 n ^« (this can be easily verfied by expanding 



10 

the expression of the total mass) . This analogy provides an interesting interpretation 
of the chemical potential (i. In fact, the Boltzmann weight for a system rotating 
with an angular momentum L and average frequency u> is exp{— /3[H — wL]}, to be 
compared with the DNLS grand-canonical weight exp{— f3[H — (iA\\, thus yielding 
the identification of uj with (i. 

Before passing to thermodynamic studies, it is necessary to comment about the 
range of validity of the XY model as an approximation of the DNLS one. The condition 
( n <C 1 implies A„ <C Z 2 = a, i.e. T <C a 2 , because on average A 2 is equal to 
the temperature T. As we are exploring the range of large values of Z = ^/a, one 
can conclude that, the larger Z, the broader the temperature range where the XY 
model provides an accurate description of the DNLS equation. Before drawing this 
conclusion, it is, however, necessary to be more careful. In fact, the presence of a 
finite conductivity in the XY model can be traced back to the existence of (possibly 
infrequent) jumps of angle-differences across the sinusoidal potential barrier. In the 
context of Eq. (|T9")) , the height of this barrier is of the order of Z 2 = a which is smaller 
than the maximal acceptable energy a 2 (since a <C a 2 , for a ^> 1). More than that, 
the validity of the XY model extends to the high-temperature regime (here "high" 
means above a) characterized by frequent jumps. 

Finally, we compare Hamiltonians (|15[) and (|19p. The former one is valid in the 
low-temperature regime, while the latter applies for Z ^> 1 (and T <C a 2 ). Accordingly, 
they should reduce to one another for small T and large Z. This is indeed the case, 
because the next -to-nearest-neighbour interaction in (fT5|) is negligible for Z ^ 1 and, 
in the low-temperature limit, one can expand the cosine interaction in (| 19 j) around 
zero. 

5.1. Thermostatted chain 

Here, we study a DNLS chain in contact with external Langevin reservoirs in the 
large mass-density limit. It is convenient to introduce the chemical-potential difference 
5(i = {(1l — A*fl)/2 and its average value (i = {(II + (Ir)/2 = 2(Z 2 — 1). From Eq. ([7]), 
we find that (for the sake of simplicity we write the evolution equation for the leftmost 
side), 

iz x = -2|zi| 2 zi - z 2 - i 1 [2\ Zl \ 2 Zl + z 2 - (2Z 2 - 2)zi - 6(i Zl ] + ^fffL^i ■ (21) 
If we assume Z>1 and a small dissipation 7, namely 

7=^ - (22) 

where 7' is of order 00, we can transform the previous equation into those for the XY 
model, with suitable dissipation and fluctuation terms 

0i = 4Z 2 Ci (23) 

C^sm^-^-V^-^+^Ci , 

where we have neglected the nonlocal dissipation z 2 and another term proportional to 
(1, as they are both negligible in this limit. By turning to the above defined rescaled 
momenta A„, we obtain 

<t>i - Ai (24) 

Ai = 4Z 2 sin(0 2 - 0i ) - 7' ( Ai - 5(i) + y/4r/'T L £1 . 
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The equations for the rightmost site of the chain are the same as the above except for 
6 \x and T L that are replaced by —S/j, and T R , respectively. Notice that the chemical- 
potential-differences ztSfi applied at the boundaries of the DNLS chain correspond to 
constant torques ±7'<5/x (i.e. a sort of external pressure) for the XY model. 

Eq. ([24]) without external torque and noise has been investigated in [33] . A closer 
version of the model ([24]) has been studied by Olla et al. [34], who, however, have 
introduced the constant-force term without relating it to the presence of the external 
bath. As a result, the external forcing does not vanish in the small coupling limit, 
with some weird physical implications for the microcanonical simulations due to the 
possible occurrence of an unbounded acceleration (if 8[i is large enough). 

5.2. Fluxes 

The continuity equations for mass and energy densities allow determining explicit 
expressions for the corresponding fluxes, 

3n = % \ z n z n-\ ~ z n z n-l) (25) 

Jn = ZnZ* n _ x + Z* n Z n -l ■ (26) 

In the large-mass limit, the leading terms read 

f n = - 2Z 2 sin(0„ +1 - <j> n ) (27) 

jl\ = -2wZ 2 sin(0 n+ i-0„)-4Z 2 n sin(0 n+ i -</>„) , (28) 

where the variables ( n have been expressed in terms of 4> n . Notice that the simple 
symmetric form of the second equation has been obtained by adding to (|26l) the 
quantity — Z 2 (<p n — <fi n -i) sin(0„ — n _i), whose average is zero in a stationary state. 
Eq. ([27]) is just the momentum flux of the XY model, i.e. its local force. The term 
proportional to <p n in Eq. (|28p has the typical structure of the energy flux in the 
XY model: it is nonzero only at finite temperatures. The first term is a coherent 
contribution that results from the fact that the oscillators rotate with an average 
common frequency to: it survives in the zero-temperaure limit, when j' r l = fij^, so 
that the heat current j^ — \ijn = 0, i.e. there is no heat transport and no entropy 
production. At low temperatures, Eqs. (|27I28|) describe the Josephson effect, where 
the chain amounts to a single junction in between two superfluids |26) . The mass 
current is proportional to the phase gradient and is independent of the system length 
N, i.e. it provides a ballistic contribution. 

6. Comparison with numerical simulations 

An important difference between the DNLS equation and oscillator chains (like the 
Fermi-Pasta-Ulam or Klein-Gordon models) is that its Hamiltonian is not the sum of 
kinetic and potential energies. Therefore, it is not obvious how to directly monitor the 
temperature T and the chemical potential fi in actual simulations. The only general 
approach we are aware of is based on non-local microcanonical expressions 35 , 10 ) 
which, at least require averaging over small subchains (see [TU] for details) . 

On the other hand, the DNLS/XY correspondence provides an effective definiton 
of temperature. In fact, from (fT^)) it is straightforward to define a kinetic temperature 
in the standard way 

Txy = (A 2 ) (29) 
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Figure 5. Comparison of the XY kinetic temperature Txy with T. Simulations 
are performed using Langevin heat baths. The dashed line correspond to 
Txy — 2T, which should hold in the limit of large fi, where the XY approximation 
is valid. The heat baths are coupled at the boundaries of a DNLS chain with 
N = 50. Txy an <i T are measured on a subchain of 30 lattice sites. 



One can easily check that in the limit of large Z, the kinetic temperature Th, 
defined in (|16p . reduces to (|29|) . Furthermore, if we compare the stochastic term in 
Eq. (|24j) with the one imposed by the fluctuation-dissipation theorem and commonly 
used in the Langevin equation for oscillator models, ■\/2 r yT (see [HE]), we can conclude 
that our definitions imply Txy = 2T (the factor 2 is just a consequence of the choice of 
the transformation of variables, and has the same origin as the proportionality factor 
J in (QU ) . 

In Figure[5]we compare T, defined as in |10j . with Txy for an equilibrium setting, 
i.e. external reservoirs at equal temperature. Txy is computed by evaluating, in the 
same simulation, the average of the („ defined in (p~0|) . The data clearly show that, by 
increasing the chemical potential /x (i.e., by increasing Z), the range of values in which 
the two temperatures coincide increases, as expected from the previous considerations. 

Finally, we have tested the validity of the XY approximation in a non-equilibrium 
stationary regime, by imposing a temperature difference AT and a small chemical 
potential difference Sfj, with respect to a reference equilibrium state (T, /x), where the 
XY approximation holds. As expected, Fig. [5] exhibits a nice agreement between 
the two definitions of the temperature profiles. A further test of the validity of XY 
model is presented in Table Q] where the full DNLS fluxes are compared with the ones 
reconstructed through the XY approximation (see Eqs. ([27| ). 



7. Conclusions 



In this paper, we have introduced suitable stochastic heat baths which can control 
both the temperature and the chemical potential in a DNLS model. At variance 
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T, T 




Figure 6. Comparison of the XY kinetic temperature profile Txy (red dashed 
line) with T [10^ (black solid line) in a nonequilibrium steady state. Simulations 
are performed using Langevin heat baths with parameters (Tl = 2, fi^ = 9.6) 
and (T R =4, f i R = 10.4). 





DNLS 


XY 


3a 
jh 


-0.234 
-1.40 


-0.208 
-1.28 



Table 1. Comparison of the exact DNLS fluxes (first column) and the ones 
reconstructed by means of eqs.l l27l) (second column) for the nonequilbrium profile 
described in the caption of Fig. [6] 



with the Monte Carlo dynamics [TO] , the Langevin equation allows establishing direct 
relationships with standard chains of oscillators. In the low-temperature limit, we 
obtain a quadratic approximation which allows describing the transport through the 
thcrmalization of linear excitations (the Bogoliubov modes), in the same spirit as 
harmonic oscillators in interaction with Langevin baths |27j . In the large mass-density 
limit, the relevant contribution to the DNLS evolution is due to the phase dynamics, 
which is described by a suitable XY model in contact with standard Langevin heat 
baths and an additional external torque. The corresponding energy flux turns out to 
be the sum of two different contributions, one due to the phase gradient associated 
with the torque, the other due to angular-velocity fluctuations. This indicates that 
transport in this region has an almost ballistic component, in agreement with previous 
observations [TO] . 

The two approximations are useful also because they allow introducing simple 
kinetic expressions for the temperature, that can safely replace the microcanonical 
definition. This is of practical importance, considering that the microscopic definitions 
of T and \i are pretty much involved for a non separable Hamiltonian, like the DNLS 
one [33] , 

An obvious extension of the present work would be to consider the DNLS model 
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on two-dimensional lattices. In this case, the correspondence with the XY model would 
predict the possibility of observing the transition from normal to anomalous behavior 
of transport coefficients at the Kosterlitz-Thouless-Berezhinskii transition [351 I3Z] • 
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Appendix: Derivation of the Langevin equation 

In this appendix we derive Eq. (O by following the system-bath coupling approach [38] . 
In analogy with what done for harmonic lattices [39] , we consider a complex oscillator, 
described by the dynamical variable z, linearly coupled with a bath of independent, 
complex harmonic oscillators described by the Hamiltonian 

H B = Y,{^ a uK\ 2 +^ h u\K\ 2 + [Kz{al + b u ) + c. C .]} , (30) 

V 

where we have introduced two different species of oscillators, corresponding to the 
two sets of frequencies u>® and oj^, while K v are the bath-system coupling constants. 
In order to preserve the global symmetry of the system with respect to phase 
transformations, we impose a second conservation law, 

^ B = £(M 2 -IM 2 ) • (3i) 

V 

The function Ab is the generator of phase transformations of the bath variables. It is 
easy to verify that the transformation generated by Ab + \z\ 2 , 

a v (s) = e js a„(0) 

&„(«) = e- %s K{0) 

z{s) = e is z(0) , 

leaves the Hamiltonian Hb invariant. An example of heat bath satisfying these 
conditions is given by a complex d'Alembert equation, Ocp(x, t) = 0, for which the 
quantity Ab represents the total (conserved) charge of the field. The equations of 
motion generated by (|3H)) are 

ia v = -wja„ - K*z 

ib v = -u v b v - K v z* 

iz = f(z) - ^2 K A a v + K) > 

V 

where f(z) accounts for the deterministic part of the evolution of z, not included in 
Hb- The first two equations can be formally solved, yielding 

o„(t) = a„(0) e"^' + iK* v [ e *" ;( *-*' ) z{t') dt' 

Jo 
ft 

b„{t) = 6,(0) e"'' + iK v / e^^-^z*^') dt' . 

Jo 

By then substituting into the equation for z, we obtain 

iz = f(z) -if G(t- s)z(s)ds + F{t) , 
Jo 
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where the noise term F(t) and the dissipation mention G(t) are defined as 



F(t) = ~Y J K u Uu{0)e 



6*(0)e" 



(32) 



G{t) 






\K U 



ioj^t —iuj^t 



(33) 



By now imposing a gran-canonical equilibrium distribution P ~ exp[— /3(Hg — /j,Ab)] 
for the bath of oscillators (where f3 — 1/T is the inverse temperature) [53], we find 
that the correlation functions of F(t) read 

(F(t)F{t')) =(F*(t)F*(t')}=0 

(F(t)F*(t')) = Y,\ K -\ 2 W^'^ (K(0)\ 2 ) + e'^-^ (\bv(0)\ 2 ) 



E 



K u 



Ju>l(t-t') 



,-iu b Jt-t') 



P(u° - n) p(u) b + n) 



(34) 



where, in order to have positive definite statistical weights, we have also to assume 
wj > fi and ur v > — /i. In the thermodynamic limit the sums over the index v in (J33]) 
can be replaced by integrals. Accordingly, we can rewrite Eq. (1331) in the form 



G{t) 



-f-oo 



dtuG a {u)e i 



dwG b {uj)e' 



(35) 



where G a ' (u>) = p a ' (lj)\K(ui)\ are two positive definite functions and p a ' (w) the 
corresponding density of states that we assume to be smooth functions. By following 
the same approach, Eq. (I34|) writes 






« w ^7 ^ 

- M P( w + W 



(36) 



which is a kind of fluctuation-dissipation theorem [39] where the Bose-Einstein 
distribution has been replaced by the Rayleigh- Jeans one. 

The corresponding generalized Langevin equation is not very practical, since it 
is non Markovian. We have nevertheless the freedom to choose the coupling and 
the density of states of the bath. The spectral properties of the process F strongly 
depend on the behaviour of G close to the ground state and may also display long-range 
correlations. To understand this point, consider the example in which G a ' b (uj) = 7. 
This choice yields a spectral density of F(t) which is logarithmically divergent close 
to the ground state frequency, thus defining a non-stationary process. The simplest, 
nonsingular case is obtained by choosing 



G» = ^(. 



M ), G\u) = ^{lo + h) 



2tt 



In this case F(i) becomes a complex white noise 



(F(t)F*(t')) = l6(t-t?) 
while the dissipation function is 

f + OC 



G(t) = 



1_ 
2tt 



duo (w - fi)e lu;t = -7 



i-S(t) + fi5(t) 



(37) 
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The full dissipation term is therefore 

-i G(t- s)z(s)ds^-jz(t)+ijfj,z(i) , (38) 

Jo 

and the resulting Langevin equation corresponds to a noisy, driven, complex Ginzburg- 
Landau equation 

(i + 7 )i = f(z) + ij/j,z + F(t) . (39) 

In the weak coupling limit (7 <C 1), the equation can be further simplified. By 
multiplying by (1 + i-y) and neglecting terms 0(j 3 ^ 2 ), one obtains 

iz={l + ij)f(z) + ij f iz + F{t) , (40) 

which has the same structure as Eq. ([5|). 
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